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ABSTRACT 

We have developed a spectrum synthesis method for modehng the UV emis- 
sion from the accretion disk from cataclysmic variables (CVs). The disk is sepa 



Wade &: Hubeny 



rated into concentric rings, with an internal structure from the 
disk-atmosphere models. For each ring, a wind atmosphere is calculated in the 
co-moving frame with a vertical velocity structure obtained from a solution of 
the Euler equation. Using simple assumptions, regarding rotation and the wind 
streamlines, these ID models are combined into a single 2.5D model for which we 
compute synthetic spectra. We find that the resulting line and continuum behav- 
ior as a function of the orbital inclination is consistent with the observations, and 
verify that the accretion rate affects the wind temperature, leading to correspond- 
ing trends in the intensity of UV lines. In general, we also find that the primary 
mass has a strong effect on the P-Cygni absorption profiles, the synthetic emis- 
sion line profiles are strongly sensitive to the wind temperature structure, and an 
increase in the mass loss rate enhances the resonance line intensities. Synthetic 
spectra were compared with UV data for two high orbital inclination nova-like 
CVs — RW Tri and V347 Pup. We needed to include disk regions with arbitrary 
enhanced mass loss to reproduce reasonably well widths and line profiles. This 
fact and a lack of flux in some high ionization lines may be the signature of the 
presence of density enhanced regions in the wind, or alternatively, may result 
from inadequacies in some of our simplifying assumptions. 



Subject headings: accretion disk, accretion disk winds — nova-like, cataclysmic 
variables — stars: mass loss — ultraviolet: general 
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INTRODUCTION 



Cataclysmic Variables (CVs) are semi-detached binary systems in which a Roche-lobe 
filling main sequence star (secondary) transfers matter to a white dwarf (primary). In the 
case of non-magnetic systems the mass transfer is made through an accretion disk. In 
Nova-like and dwarf novae in outburst, the disk is the main source of radiation. The inner 



disk emission dominates the ultraviolet (UV) spectrum ( IWarnerl 119951 ). 



Several efforts have been made to understand the spectral features of accretion disks 



within different physica . 



Lynden-Bell fc Pringle 



frameworks, 



1974 : 



Tylenda 



rom s i mple composi te black-bodies ( iLynden-Bell 



1969 



1977 



Pringle 



198ll ) to more realistic atmosphere-disk 



emission m odels that tak e into account line blanketing, Doppler broadening and limb 



darkening ( iLa Dous 



1989 



Diaz et al. 



1996 



Wade fc Hubeny 



19981 ). The influence of the 



white dwarf (WD), disk rim and orbital ph ase on the spectrum have also been studied 



( iLinnell fc Hubeny 



1996 



Linnell et al. 



20071 ). These attempts to 



features have been widely 



1999 



Nadalin fc Sion 



t ested with UV data (e.g. 



2001 



Engle fc Sion 



2005 



Long et al. 



Linnell et al. 



reproduce UV continuu m 



2008 



994 



Diaz fc Hubeny 



and there have been 



some successes in reproducing some photospheric spectral features, for example the 
absorption proflles of Si III AA1298,1304 and Cll A1335, and the light curve through eclipse. 
However, these models have deflciencie s — they are unable to reprodu ce the flux level and 



continuum color at the same time (e.g. 



Wade 



1988 



Puebla et al. 



20071 ). Another signiflcant 



problem with these models is that even though they reproduce some photosphere spectral 
characteristics and some absorption proflles, they fail to reproduce the optical and UV 
emission lines that are observed for CVs. 



Since the pioneering 



(^Tylenda 



1981 



Williams 



studies to model the emission lines from dwarf novae in quies cence 



1980 



Williams fc Ferguson 



1982 



Williams fc Shipman 



1988|), new 



physical scenarios have been proposed to explain the strength and the orbital behavior 
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of emission lines. The P Cygni profiles of Civ AA1548,1551, Siiv AA1393,1402 and Nv 
AA1238,1242 observed in low to intermediate orbital inclination systems in high state, 
provide convincing evidence of mass loss through strongly accelerated winds (terminal 
velocities ranging from 2000 to 5000 km s~^). The observed profiles resemble UV spectral 
features observed in the O B stars, and this has led to th e idea that these outflows could 



also be line-driven winds ( ICordova fc Mason 



1982 



19851 ) driven off the accretion disk. 



However, some P Cygni profiles in CVs show differences from those observed in OB 
stars. For example, the deepest flux level in the blueshifted absorption component is 
closer to the transition wavelength, while for hot stars it is close to the wind terminal 



velocity (e.g. 



Shlosman fc Vitello 



1993 



Prinja fc Rosen 



19951 . and references therein). In 



addition, a correlat ion between disk wind signatures and the accretion rate has been found 



(IPuebla et al. 



20071 ). For systems with low accretion rate (e.g., dwarf novae in quiescence) 
the absorption component is absent or indiscernible in an asymmetric emission profile, while 
for high accretion rate systems, a deeper blueshifted absorption component is often seen. 

Studies on eclipsing systems indicate an axi-symmetric wind instead of spherical 
geometry and a highly stratified ionization structure of the emission line region in CVs 



(IKing et al. 



1983 



Drew 



1987 



Mauche et al. 



1994 



Mason et al. 



19951). The blueshifted 



absorption component gradually weakens as the orbital inclination increases. Thus, for high 
inclination systems the spectrum is mainly dominated by quasi-symmetric emission lines. 



Mason et al 



(119951 ) showed variations on the line profile and intensity through the 
eclipse in UX UMa, and concluded that the lines are less eclipsed than the continuum, 
especially Civ AA1548,1551. This suggests that the size of the emission region is comparable 
to that of the secondary star, and the line emitting region is larger t han the continuum 



emitt ing region. Similar results were found for RW Tri and DQ Her (ICordova fc Mason 



19851 ). More recently high-time resolution far-ultraviolet (FUV) observations of UX UMa 



- 5 - 



reveal strong and blueshifted Cm, Nlll and O VI absorption lines (e.g. iFroning et al.ll2003[ ) 
that suggest active mass loss. 

The first effo r t to u nderstand disk-wind lines quantitati vely was ma de by 



Cordova fc Mason! ( 1l982l ) who used pure stellar wind models (jOlsonI 119821 ) to fit the 



deep absorption components of some lines from the dwarf nova TW Vir in outburst. These 



authors adopted a spherical velocity profile centered on the whit e dwa r: 



the mass loss rate as 10 ^ to 10 ^ of the accretion rate. Later, 



Drew fc Verbunt 



Drew 



, and estimated 



(Il986[ ) (see also 



19851 ) studied the effect of an extended continuum source on line profiles. 
They pointed out the strong influence of the mass loss rate on the line strength and the 
weak effect on the blue wing absorption. These authors also found that the more dense the 
wind, the bluer the P Cygni absorption minimum. Finally, they applied their theory to lUE 
observations of RW Tri and UX UMa, and calculated a mass loss rate >10"^° M0 yr~^. 



A second generation of models was developed by iShlosman fc Vitelld (jl993[ ) (hereafter 
SV93). To generate proflles for Civ, they developed a kinematical model with a bi-conical 
geometry, and with a stellar velocity law along tilted stream lines. A n extended blackbody 



continuum source was constructed from a 



Shakura fc SunyaevI (119731 ) flat disk, the white 



dwarf, and the boundary layer (BL). These authors assumed an isothermal wind with an 
ionization structure calculated through the UV photoionization-dominated media (H, He 
and CNO) and used the Sobolev approximation for radiative transfer. They found a series 
of proflles that show the rotation eff ects on line shape, and the strong influence of wind 



geometry ( IVitello fc Shlosman 



19931 ). They estimated the mass loss rate to be within ~4% 



to 15% of the accretion rate. Later, they showed the i mpact of wind rotation on the proflle 



for high orbital inclination systems such as V347 Pup ( IShlosman et al. 



19961 ). 



Another effort to understand the wind pr operties through lin e synthesis was made 



using Monte Carlo radiative transfer methods. 



Knigge et al 



( I1995I ) assumed a stellar wind 
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velocity law, an isothermal gas, and that the line is formed by pure scattering. These 
methods helped to analyze the profile structure for a single line in more detail. After 
that, Monte Carlo methods were also used to calculate the ion and temperature structure 



of a bi-polar wind using the Sobolev approximation. 



the fu 



1 Cy spectrum in a broader wavelength range ( iLong fc Knigge 



'his was the first effort to model 



2OO2I). However, 



Long] (120061 ) pointed out the difficulties in reproducing more than one observed line wit h 



the same wind parameters. More recently, this method was improved (jSim et al. 



20051) 



allowing t he formation of recom bination lines and applied to RW Tri and UX UMa HST 



UV data f lNoebauer et al 



20101) • These models yield interesting and encouraging results 
about line profiles and disk wind structure and their results are compared with ours bellow. 
Particularly, previous models do not self-consistently compute the disk-wind interface. In 
this work we attempt to do it. 

The aim of this work is to better understand the ionization structure, geometry and 
physical parameters of the line emitting region in CVs. To achieve this goal we have used 
disk- wind models, that include detailed wind physics ([|2]), to simultaneously model several 
line features observed in the UV spectra of CVs. The dependence of synthetic spectra on 
parameters such as the white dwarf mass (Mwd), accretion rate (Mq), mass-loss rate (M^), 
orbital inclination (i), and wind geometry was investigated. A description of the method, 
including the calculation of the velocity laws, wind structure and spectrum synthesis is 
presented in ^ Extensive simulations are shown in §31 A comparison of synthetic spectra 
with UV and FUV data for the high inclination CVs RW Tri and V347 Pup is made in §11 
A general discussion of our results is given in §51 while a brief summary and conclusions are 
presented in §61 
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2. THE METHOD 

2.1. The Disk- Wind Model 

The observed UV hne profiles of many high mass transfer CVs resemble those observed 
in O stars winds. Due to the similar spectral characteristics, it is feasible to suppose that in 
CVs the emission lines are formed in a wind ejected by the system and, similarly to stars, 
this outflow is driven by line radiation. 

Actually, the disk wind is a complex 3D structure, but as was pointed out in section [H 
the wind seems to have an axi-symmetric geometry. From this fact, the models should use 
at least an 2.5D (2 spatial directions plus rotation) approximation, similarly 2.5D radiative 
models need to be used to predict the observed disk-wind spectra. However, a complete 
2.5D hydrodynamic plus radiative transfer calculation would be much more complicated 
and time consuming than the ID modeling. Since many basics of CV wind structure and 
acceleration are still not understood we have adopted an intermediate approach. We use 
ID models to compute the atmospheric structure (including the deep layers of the disk, 
disk photosphere and wind) - that is, the temperature structure, the ionization structure, 
and the level populations in vertical sense. We then combine a series of ID models to 
create a 2.5D model, and use this 2.5D model to predict synthetic spectra. As we will show 
below, in this approximation the temperature in the extended wind region is controlled only 
by the radiation within the ID model corresponding to the r radius. This would not be 
desired, taking into account that the wind temperature far from the disk is controlled by 
the radiation from the whole disk. Despite these simplifying assumptions we got reasonable 
good UV line profiles when compared with observations (section H]). One advantage of our 
approach is that we correctly treat the photosphere-wind transition region as is described 
in section [2.2.11 
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The model disk atmospheres and synthetic spectra are calculated using an arrangement 
of concentric wind atmospheres. The base of each atmosphere is located at the disk plane 
at a distance "r" from its axis of symmetry. Each wind atmosphere has its own internal 
disk structure, disk photosphere and vertical wind. The atmospheres are placed side 
by side, thus forming the 2D structure of disk, photosphere and wind. The atmosphere 
structure calculus requires a vertical velocity field, which is obtained consistently with the 
physical conditions in the disk, such as gravity, radiation and density for a wind within 



Castor et al. 



fll975l ) (hereafter CAK) framework, as explained in section I2.1.1[ Once the 
atmosphere structures are calculated the opacity and emissivity in the co-moving frame 
(CMF) are computed for each depth in the vertical grid. These values are then stored on 
a two-dimensional (and finer r, z) 2D grid through interpolation. This grid covers from an 
inner radius rj, to an outer radius Vf as it is shown in figured^. The vertical size of the 
wind is set equal to Rdisk- The region between rj and Vf embraces the radii where the wind 
is expelled. Each of wind atmosphere is calculated in the vertical sense using the method 
described in section 12.2.11 that is widely used in the of study hot star winds. We obtain 
thus, a cylindrical wind with a vertical velocity field. 

It is necessary to introduce two more components to velocity field, a radial component 
Vr and azimuthal component V^, that together with will bear a 3D velocity field. We 
arbitrarily set the Vr component in order that the projection streamline describes an 
hyperbolic trajectory with an asymptotic angle /3 (see figure [2]). The streamlines don't 
follow ID vertical models, their bases are on a point of the disk plane with radius To 
corresponding to an hyperbolic trajectory that passes through the (r, z) point with and 
Vr- The Vtf) component was calculated in order to conserve the specific angular momentum 
from the wind base (keplerian) and to simulate the divergence of the wind streamlines. The 
base stream radius Tq was calculated from the trajectory equation for each (r, z) point: 

z^ = tan''(3{r^ -rl) (1) 
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Thus the velocity field is given by 

-—^7-^Vz(r, z), if z < r tan(/3), 
Vr{r,z) = { "-'^-"^^^ (2) 

4^V;(r,.) ifz>rtan(/3) 
and the specific angular momentum conservation (ro a/GMwd = rVfj,{r, z)) leads to: 



V^{r,z) 




r-2 tan2(/3) 



1/2 



^ if ^ < rtan(/3), 
if 2; > r tan(/3) 



(3) 



where K is the radial velocity, is the azimuthal velocity, /3 is the asymptotic aperture 
angle of the streamlines (figure [2]) and Vz is the velocity law calculated as shown in §2.1.1[ 
The parameter h is the corresponding value of z of the last vertical grid point such that 



V^>0. This is to avoid discontinuity beyond z=r ta.n{P). In t 
Such a value corres p ond t o an average of the values found in 



lis work we use =45°. 



Proga et al. 



fll999l ) and 



Pereyra fc Kallmanl ( l2003l ). It was kept fixed as we needed to constrain the number of free 



parameters to explore. 

In the spectral synthesis we have introduced additional parameters to modify this 
rigid geometry in order to explore their infiuence on the spectral features. The figure 
[T}3 shows those parameters and their meaning. The aperture angles 61 and 62 limit the 
regions which will be taken into account in the final radiative transfer calculus. There is 
no matter within the cone defined by 61 and outside the cone defined by 62- In the last 
case, the region inside the cone defined by 62 and beyond is filled by the same structure 
corresponding to the last wind atmosphere calculated on r/. The Tc radius also limits the 
regions that eject wind, but keeping the wind structures considered above the line defined 
by 02. There is not a consistent reason for these approximations, we introduce them in 
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order to ra ake our analysis more flex i ble and to mimic the well studied geometry of disk 
winds (e.g. 



Shlosman fc Vitello 



1993 



Long fc Knigge 



20021 ). Using those 



were able to follow the geometry predicted in hydrodynamic models (e.g 



Pereyra et al. 



2000 



Pereyra fc Kallman 



parameters we 



20031). 



Proga et al 



1998 



It is clear that for regions close to the symmetry axis the corresponding streamlines 
wouldn't hit the disk. Those streamlines will be empty (or fllled) depending on the 6i value, 
this is parametrized to simulate the kinematics of the wind in these regions without solving 
the wind dynamics in 2.5D. Also, from this approximation the mass is not exactly conserved 
trough the wind streamlines, specially in the regions far from the disk photosphere. Despite 
of this, we can calculate an approximate value for the total mass-loss rate, M^^, by: 



= Ait m{r)rdr (4) 

where m(r) is the mass loss flux from disk that is calculated from Euler equation as it is 
shown below (section 12.1. ip . An additional factor of 2 has been included to account for 
both sides of the disk. In our approximation the calculated mass loss value will be slightly 
different than the real value that would produce the line emission. This due to the influence 
of ^1, and rc on the volume size that is taken into account in the spectral synthesis. 
The more these parameters are away from a cylindrical model, the farther this value of 
mass loss will be off. However, an increase in 9i is roughly compensated by an increase in 
^2- Considering the range of parameters used in this work, we can estimate that in worst 
cases the departure can reach ~20 %. 

Using the resulting emissivity and opacity values, mapped into the observer's frame, 
the radiative equation is solved exactly along rays that are deflned by a set of impact 
parameters distributed concentrically in the disk plane. By integrating outgoing speciflc 
intensities over the solid angle comprised of the disk plus wind we obtain the synthetic 
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spectrum. The details of each step are described in the following sections. 



2.1.1. The vertical velocity laws 



The velocity law is required for defining the density profile and radiative transfer 
in the disk wind. It is calculated by solving the Euler equation for a vertical wind with 
the conditions found at each ring location in the accretion disk. It is important to note 
that unlike stellar atmospheres the gravity at a given radius increases with height inside 
the atmosphere up to a maximum at ^^0.7. Gravity then decreases at greater heights, 
following a ^ law at large distances above or below the disk. The flow base is at the 
disk plane and extends out to a distance equal to the disk rad i us. H ere, we solved the 



momentum equation in the context presented by 



Pereyra et al 



)3) 



(120041 ). This approximation 



has already been used to study hydrodynamic properties o: 



[e.g. 



Vitello fc Shlosman 



19881 



Pereyra et al. 



1997 



an isothermal wind can be expressed as follows (jPereyra 



Perevra 



disk winds in previous works 



20051): 



20051 ). The Euler equation for 



s\ du / 1 du\'^ 2s da 

1 1 -T = -9^ + fa[—— ] 

00/ dx \m dx I a dx 



(5) 



and following the 



Pereyra et al. 



( 20041 ) and iPereyral ( 20051 ) reasoning for of a vertical 



disk- wind we found that: 



-VI- 



VO 



m 
rho 

/ 
9 



z 
,2 



V. 



o 



X = arctan(2/r) 



sec^(x) 
0,2 



rh{r)/mo 

a (1 — ay " %' 



{l-a)/a 



^faT^{r)k ( 1 



■re/ 



ll — o ' ^ ' 



sm x cos X — 



— cos''(x) 



cGM 



WD 



(6) 



In the above, the effective temperature used in this worlc follows the standard accretion disk 
model radial profile: 



T{r\ 



3GMwdM, 



1/4 



WD 



R 



WD 



3/4 



R 



WD 



1/2- 



1/4 



(7) 



see 



Shakura fc Sunyaev 



1973 



Lynden-Bell fc Pringle 



I974J ) where r is the distance from 



the disk center to respective annulus measured at disk plane, Rwd is the compact star 
radius, o is the Stefan-Boltzmann constant, G is the gravitational constant, is the 
mass accretion rate and M wd is the WD mass. Assuming a K pure carbon WD of 



-0 7 

Hamada fc Salpeterl (119611 ). Rwd oc Mwd ' and the maximum disk temperature T 



scales as ~Mwd • Furthermore, v is the velocity along a streamline that begins on 
the disk, a is the geometrical factor that takes into account the streamline divergence 
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f lVitello fc Shlosmanlll988[ ). Wo=V2GMwD/r is the escape velocity at the base of the wind, 



b=^JkBT{r) I iivfiii is the sound speed, is the Boltzmann constant and /tr'^-^=Q.325 cm' 



gr ^ is a reference value for the electron scattering opacity ( iHowarth fc Lamers 



19991). We 



have used an atomic weight /i=1.28, similar to the solar value, and mn is the proton mass. 
The CAK parameters a and k are physically related with the line opacity distribution and 
with the ratio of optically thick lines to optically thin lines considered in the line force. 
Thus for a=\ only thick lines are taken into account, whereas for a=0 only thin li nes are 



consid ered. These parameters are not independent of each other a s pointed out 



Gayleyl fll995l )) show a 



by 



Gayley 



fll995l ) (see also CAK). The measured a and k values (Table 1 in 
power-law distribution of lines and power law relation between these parameters. In this 
work we use this distribution for a range of a between 0.9 and 0.5. The function g evaluates 
the gravity force corrected for electron scattering acceleration and normalized by GMwd/''"^- 



The mass loss flux rh(r) is normalized by rho, which in turn de pends only on 



characteristics. The parameter iho is the equivalent to McAir in 



Pereyra et al. 



ocal disk 



mm . 



The function (^[x) represents the disk contribution to the radiation fleld, at a particular 
wind point. This radiation comes mainly from the disk surface and is used to calculate 
the line radiation force. Such a this force is normalized to the electron scattering force 
th rough the 7n p a ramet er. In this work, we mainly used the analytical relation proposed 



by 



Pereyra et al. 



(119971 ) that takes into account the ra d iation contribution from the whole 



disk. Also, for some cases, we used the 



Pereyra et al. 



(l2004f ) rela tion that prescribes / 



Pereyra et al. 



mm 



models for inner regions and O models for outer ones. Following 
the / models try to mimic the vertical flux distribution of the inner standard disk model, 
and the O models try to mimic the flux distribution of the outer regions. The former falls 
as [1 — (^)^]^^ because most of radiation is coming from region immediately below. In 
the O approximation the flux flrst grows with height as the wind see the inner (hotter) 
regions, then it decreases in the same sense as in the / models. The match criterion for 
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both models is the same terminal velo city. We designate Ci (eq. [8]) for 



parameterization and (2 ioT the one of 



equation [S] behavior is presented by 



Perevra et al. 



Pereyra et al. 



Pereyra et al 



fll997h 



mM) ( eg. ED. A detai led analysis of 



fl2004h and 



Perevral (120051 ) . 



Ciix) = sin"(x) cos^(x) 



(8) 



C2(X) 



4+tan2(x) 



/ case 



l+2tan(a;) 2 f \ 

-TTT — rvcosix) Cease 

l+tan(x) ^ ' 



(9) 



In this work we iteratively solve the Euler equation to find the critical point and the 
mass loss corresponding to the critical solution (CAK). The boundary conditions fo r each 



iteration were set by calculating an atmosphere disk model (jWade fc Hubeny 



1998h and 



a velocity law that follows its density structure. This velocity law is computed from the 
disk plane to an outer optical depth rij=10~^ (where is the Rosseland optical depth). 
Then we find the height Zc where t^{^Z(^=\ and the corresponding vertical velocity and 
acceleration v'^ at this point. Equation [5] is solved with the boundary conditions vi^z^^j = 



and 



dv I 
dz ' 



v'^. The theory of stellar atmosphere winds shows that the Euler equation has 



two solutions, the one of them, the lower branch, goes from subsonic region to the Parker 
point (in our case where h = ga — = 0). The upper branch solution goes from infinity 
to the sonic point (where Vz = b). Thus the re gion where the two solutio ns co-exist is found 



I999I ). The only solution 



between the sonic point and the Parker point ( iHowarth fc LamersI . 
that can start subsonic and grow smoothly through the sonic point to infinity is the critical 
solution that connects the lower branch solution to the upper branch solution at the critical 
point (CAK). The lower branch solution, even supersonic, fails when it reaches the Parker 
point, but in our simulations this point has never been reached. That is the reason why we 
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adopted the lower branch solution in our spectrum synthesis. These solutions bear lower 
terminal velocities that are more comp atible wi t h the observations. The critical solutions 
are found using the criteria exposed by iPereyral (120051 ) for the existence and continuity. 



The dependence of the terminal velocities on a, as well as the dependence of the mass 
loss flux, are shown in Figure [31 These values were calculated using the "C2" function with 
MwD= 0.8 and M©, Ma=lxlO~^ Mq yr~\ It was found that the critical point position is 
strongly dependent on a, as are the terminal velocities. The mass loss flux also is dependent 
on a and /c, as evident from Equations [5] and El it has also a close relation with the effective 
temperature at the base of each wind model. This relation is similar to that observed in 
hot stars, as well as the relation between the total mass loss rate and luminosity. 

For our disk models, the mass-loss rates (integrating up to r=0.3 R0=2xlO^*^ cm) 
are 6.0x10"^° M0 yr"-*^, 5.0 xlO"^-*^ Mq yr~^ and 3.4 xlO"-*^^ Mq yr"-*^ for decreasing a 
values from 0.9 to 0.5. These mass-loss rates are consisten t with 2D and 3D hydrodynamic 



simulations (jPereyra et al. 



1997 



Pereyra fc Kallman 



2OO3I ) and represent ~ 0.04 % to - 10 



of the accretion rate. 



2.2. Spectral Synthesis 

Once the wind model grid is calculated we are able to compute the CMF emissivity 
[jiy) and opacity (x;^) at each point. The opacities and source functions (S,^) in the CMF 
are interpolated into a finer grid in the zr plane. Additionally, they are mapped into the 
observer's frame. In order to synthesize the spectrum, a new grid of log-spaced impact 
parameters on the plane of the disk is set from pmm=l-05RwD to 'Pmax = -R(iisfc[tan(z) + 1], 
where i is the orbital inclination. In addition, a linearly spaced grid of azimuthal angles 
is set. A ray goes to the observer from each point of the 2D grid thus formed. Radiative 
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transport through the inner disk and wind is calculated using points along each ray that 
are spaced 1 km s~^ in projected velocity. For impact parameters greater than tq and 



lower than Kdisk the corresponding value o. 
atmosphere at this r adius as describe d by 



l,Az = 0) is equal to t he emissivity of a disk 



the limb darkening (jPiaz et al 



Wade fc Hubenyl (119981 ). taking into account 
19961 ). Thus, for impact parameters pmm < p < rc 



calculate the radiative transfer consistently from disk plane through the photosphere and 



wind. For impact 



Wade fc Hubeny 



param eters Vc < p < Rdisk the disk is treated as a surface that emits as a 



(jl998[ ) disk-atmosphere. Finally, for impact parameters Rdisk < P < Pn 



l,{z = 0)=0. 

The radiative transfer is calculated using a broad frequency grid in CMF and a 
frequency grid in the observer frame. For each point in the latter grid the radiation is 
calculated at each point along the ray. This is done through the corresponding observed 
frequency in CMF, using the shift caused by the velocity field projected in the observer 
direction at each point. The projected velocity is given by 



Vp = Vr COS COS i + sin cos i + VzSini . 



(10) 



The resultant specific intensities 1^, in the observer frame, that leave the system from 
each impact parameter are then interpolated into a finer frequency grid and also in a finer 
azimuthal grid. Those specific intensities are integrated over the solid angle of the disk plus 
wind to get the synthesized spectru m of the whole system (disk+wind). Our 2D code was 



adapted from an earlier 2D code by 



Busche fc Hilherl (120051 ) . used to compute spectra from 



rapidly rotating stars that have an expanding atmosphere. The main change involved a 
switch to polar coordinates and disk-like geometry, rather than the spherical coordinates 
used in the original code. 
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2.2.1. The ID Disk-Wmd Atmosphere through CMFGEN 



code 



he method used to calculate the ID disk-wind structures is based on the CMFGEN 
Several details this method are widely described by 



Hillier 



1990 



Hillier fc Lanz 



2001 



Hillier 



Hillier fc Milled mm (see also 



20031 ). CMFGEN is used to solve simultaneously 



the statistical equilibrium equations (non-local thermal equilibrium or NLTE), the radiative 
transport equation in the CMF, and the radiative equilibrium equation. For simplicity 
(and feasibility) we assume that the perturbed radiation field depends only on the local 
populations (diagonal approach), or local plus neighboring depths (tridiagonal approach). 
The latter has generally better convergence, but requires twice the memory. The only 
difference between these methods is the model convergence. Once the convergence is 
achieved, the results are the identical. This process is repeated until convergence of the 
populations is obtained. CMFGEN allows for line coupling, level dissolution, dielectric 
transitions and electron scattering. In addition, charge exchange. Auger ionization, and the 
influence of X-rays can be treated — these processes have not generally been included in 



previous accretion disk wind models (IHillier fc Miller 



19981 ). although their influence on the 



line proflles has not been quantifled yet. 

The energy balance constraint, as expressed by the radiative equilibrium equation, is 
used for deriving the vertical temperature structure. Line blanketing plays a crucial role in 
determining the atmosphere structure and spectra and it is included. In the case of hot 
stars, line blanketing has a strong influence for A < 1300 A. A similar statement applies for 
accretion disks, especially in the inner regions which have similar effective temperatures. 
Originally CMFGEN was written to treat spherically symmetric expanding outflows; this is 
not the case in accretion disk, whose geometry tends to be bipolar. CMFGEN was extended 
to allow the calculation of wind atmospheres with a plane-parallel geometry, as needed for 



^http:/ /kookaburra. phyast.pitt.edu/hillier/web/CMFGEN. htm 
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thin disk models. In this work we use this geometry for each cylinder (disk plus wind). 



Even with a plane parallel geometry, CMFGEN still uses stellar-like parameters, 
namely, the luminosity (in L0), core radius (Rp in 10^'^ cm), external radius (in Rp units), 
mass loss rate (M^,), besides the wind velocity law, atomic models, ionization states present 
in the gas, atomic species and their abundances. In the disk case we have different local 
physical conditions at each radius. These physical conditions are characterized by the 
effective temperature, mass loss rate and the height to which the wind extends. In order to 
get the correct input for CMFGEN a fiducial core radius Rp=10.051 is set for all models, 
while the external radius is set to the core radius plus disk radius, which is set equal to 0.8 



Rij, where R^^ is the distance from the disk center to the inner Lagra ngian point 



[T^). This value is close to the truncation limit imposed by tidal forces ( lOsaki et al. 



'figure 



I993I). 



Mass loss rates are obtained from the local mass loss fiuxes, and luminosities are obtained 
from the stellar radius Rp and the effective temperature. 

The wind atmosphere models are calculated tak ir ig int o account the internal d ensity 



structure of disk as is described by 



Wade fc Hubenvl fll998h (see also 



Hubeny 



1990). In 



these models the radiation field and physical parameters deep inside disk-atmospheres 
depend on the highly uncertain behav ior of viscous dis sipation, which is parameterized as a 



power law of the column density (see 



Diaz et al 



19961 ). Once this structure is calculated. 



the corresponding vertical velocity field up to tii=1 is computed using continuity equation 
for a thin disk, p{z)v{z) = m{r). This velocity field is matched with the Euler solution as 
was explained above. The structure calculus with CMFGEN thus, goes consistently from 
disk plane passing through the disk photosphere, wind-photosphere interface out to wind 
with high z as high as Rdisk- 

The mass loss rate per unit area at the photosphere (m(r)) is calculated at each 
radius as explained in the section 12.1.11 The wind atmospheres are not calculated out 
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to the external disk radius. Instead, the maximum radius used to calculate the wind 
atmosphere depends on the effective temperature, and generally stops around ~14000 



K. The contributions to UV spectra and the inffuence on the main wind lines 



these regions are small. Accretion disk photospheric emission (jWade fc Hubeny 



jeyon d 



19981) is 



calculated for the outer regions and transfered through the inner when the rays shock it. 



3. MODEL RESULTS 

In this section we analyze a set of simulations with the aim of studying the dependence 
of the UV spectrum on the model parameters. To do this, we calculated the disk structure 
and wind structure, as explained in the previous section, in a limited parameter space 
which covers the observed one, as found in literature. A range of CAK parameters {a,k) 
that bracket the values for hot stars is explored. For a given accretion disk model, a and k 
define the mass-loss rate and the wind velocity law for each disk ring. Tables [T] and [2] show 
the parameters used in this work, with the physical parameters that characterize the wind 
listed in Tabled It is important to note that the physical parameters are not independent, 
actually only the Mwd, M^, a, k and ({x) are. The remaining are dependent, namely, T* 
depend on Mwd and M^, as well as and Voo. In this work we have fixed the /3 value and 
the ({x) function in order to not increase the number of free parameters. A brief discussion 
in the influence of the ({x) function on the "photospheric" line Cll A1335 is presented 
bellow (see section [5]). 

The primary mass (Mwd) and the mass accretion rate (M^) determine the disk 
properties, which in turn act upon the wind temperature and wind ionization structure. 
The mass loss flux is not an independent parameter - rather it is determined by the 
parameters describing the accretion disk, and the CAK parameters {a,k) through Equation 
(|5]). We also varied (see Table |2]) the angles 6i and 6*2 that limit/define the calculated wind 
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region that is taken into account for the spectrum synthesis. 62 mainly determines the 
extended atmosphere size at the coolest wind atmosphere calculated (Fig. [T]). The other 
geometrical parameters are the orbital inclination (z), the wind height (Z^ax) and the disk 
radius (Rdisk) that set the wind volume and the part of the wind hidden by the disk itself. 

3.1. Atmosphere structures 

In this section we analyze the atmospheric structure in the vertical direction. The main 
factors that determine the ionization structure and the level populations are the electron 
temperature and density. As the structures of each disk-wind atmospheres are similar to 
each other, we chose one of the calculated models (model "e" ) in Table [1] for the analysis 
— differences between models will be highlighted, as necessary. Model "e" has a white 
dwarf mass = 0.8 Mq and an accretion rate of Ma=1.0xlO^^ Mq yr^-*^ that yields 

a disk-temperature compatible with the most luminous CVs. The total mass loss-rate is 
9.3 xlO"^"*^ Mq yr"-*^, which is 1% of the mass accretion rate. While the mass fluxes vary 
by over two orders of magnitude, the outer regions of the disk (because of the weighting 
by area) make a very important contribution to the total mass-loss rate — for model "e" 
each shell (i.e., the area between ring i, i+1) contributes roughly an equal fraction [within 
a factor of 2] of the total mass-loss rate. The figure H] shows this characteristic in a clearer 
manner for model "e" . While the flux mass falls with increasing radius, the total mass loss 
(quantified by A7rm{r)rdr) tends to stabilize at large radii. 

The physical parameters of each atmosphere are shown in Table [31 In this table the 
radius, in WD radius units, is shown for each ring, as well as the local mass loss flux, the 
vertical terminal velocity and the effective temperature following equation ([7]). This table 
also shows the star-like parameters used as input, namely: the stellar luminosity (L^.) and 
the total mass- loss rate (Mj,). 



- 21 - 



The last two atmospheres, corresponding to the largest and coldest rings, were not 
calculated. There are two reasons for this. The first one is physical — the contributions 
of these regions to the UV emission are low due to their low temperature. Further, 
these regions would only contribute to the line profiles if the wind density above the 
corresponding rings is high. The second reason is numerical — we found that, using 
plane-parallel geometry, the algorithms face convergence problems in cold, low log(5') and 
low density atmospheres. 

Detailed level populations in NLTE were considered for the H, He, C, N, O, Ne, Si and 
Fe. Their ionization states are shown in Table H] with the level and super-level numbers 
used for each ion. 



3.1.1. Temperature 



In Figure E] the vertical electron temperature profiles are shown for each ring of the 
disk plus wind system. In the left panel the temperature structure for the photosphere-wind 
transition region is shown. Every temperature structure for each ring displays the 
same morphology. From deep in the atmosphere the electron temperature decreases 
monotonically towards the photosphere, where the radiation begins to uncouple from 
matter. The temperature then rises because of an indirect effect of Lyman and Balmer lines 
on the Lyman-continuum heating until a value lower but close to the corresponding effective 
temperature is a t taine d for each ring. This is a classical NLTE effect first discovered by 



Auer fc Mihalad ( 119691 ). Note that the minimum temperature happens at a different height 
for each ring, as expected from the different height of tii=1 for each disk photosphere 
structure. This is caused by the dependence of the "z" component of the WD gravity (which 



sets the atmospheric scale 



flaring (IFrank et al. 



2002 



leight in the disk) on z and r that generates the well known disk 
, pag. 87). The right panel in Figure [5] shows the temperature 
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behavior for extended wind regions for each atmosphere modeL There the temperatures are 
almost constant with height, with a value close to the corresponding T(,ff{r). This behavior 
can be explained by the lack of radiative dilution in the plane-parallel approximation. Thus, 
the average outer electron temperature is Te„~0.85Te//(_R). Actually the ratio Te^/Te//(i?) 
is smaller for inner radii (~0.75), and larger for the outer zones where this ratio is almost 
equal to 1. 



3.1.2. Ionization Structure 

Due to the isothermal nature of the outer wind, and because there is very little 
dilution of the radiation field in the plane-parallel approximation, the ionization state for 
each atmosphere model is almost constant in the wind above the photosphere. However, 
strong variations in the ionization states are found in the transition region between the 
photosphere and wind. These variations are stronger in cooler atmospheres. The ionization 
state for carbon, for ring 6 of model "e" which has an effective temperature of 30300 K, is 
shown in Figure [61 Here, the left panel shows the photosphere-wind transition zone, and 
the right panel shows the extended wind region. The former indicates a highly structured 
ionization state, in particular above r/j=l. The dominant ion changes from C"'"'^ in the 
inner photosphere to C"'"'^, and from here to C^+. This happens in a small region with 
the electron temperature attaining its minimum at r/j^SxlO"'^. Immediately above this 
point, temperature begins to rise and the C"*"^ concentration rises in the same way. When 
r/j~1.8xl0~'^ and the electron density Ne~3.7xl0^^ cm""^, it becomes the dominant ion for 
the carbon. For greater heights the ionization structure tends to stabilize. That trend in 
the vertical direction is the same for all models. 

As a consequence of the changes in effective temperature with radius, there are large 
radial changes in ionization. For instance, at the inner radii, the dominant ion for carbon is 
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the C"*"^ until the wind temperature is ^39000 K at r = 4.6 Rwd from the disk axis. For 
temperatures between 39000 K and 24000 K and 4.6 Rwd < r < 9.7 Rwd the dominant 
ion is C^'^. Beyond 9.7 Rwd? and for temperatures lower than 24000 K, the dominant 
ion is C''""*" (out to r = 29.6 Rwd)- Furthermore, if a high value for 62 is used, the last 
model atmosphere is larger, and its ionization state becomes dominant in the wind. This 
behavior is followed by the other ions in the wind, depending on the radiation temperature 
that is coming from the model disk photosphere. In the same way, the level populations 
have strong variations in the photosphere-wind transition region, and a stable behavior in 
the wind zone. Beyond tji=1, the level populations depart significantly from LTE. The 
departure coefficients increase with height up to a maximum, decreasing as gas enters the 
accelerated zone. 



3.2. Synthetic spectra 

The effects of changing the model physical and geometrical parameters on the wind 
line profiles and continuum are described in this section. We have chosen models from 
Table [T] that differ in one physical or geometrical (Table parameter to guide the analysis 
of the UV observations. 



3.2.1. Line profile dependence on physical parameters 

Models "a", "fe", "c" and "ci" have been chosen in order to investigate the influence 
of the accretion rate on the line proflles. These models have the same white dwarf mass, 
but different accretions rates. Due to the dependence of the mass-loss rate on the effective 
temperature through Equation ([5]), these models have different M^„. Consequently we have 
two effects, the change in wind density and the change in wind and disk temperature, both 
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arising from changes in the accretion rate (Eq. [7]). 

In figure [71 the variation in fine profiles with accretion rate is shown for Cm A1175 
and Civ AA1548,1551 UV fines. The accretion rate acts in a different way on each fine but 
generafiy tfie blue absorption component is deeper wfien tfie accretion rate is fiigfier. Wfien 
tfie accretion rate is fiigfi tfie mass-loss rate is also fiigfier and tfiere are more absorbers in 
tfie wind. Furtfier, an increase in tfie accretion rate increases tfie disk temperature and 
fience tfie wind ionization, and tfie relative intensities between lines cfianges accordingly. 
For instance, tfie intensity of a low ionization line like Cm A1175 increases wfien Mq 
decreases, wfiereas tfie intensity of Civ AA1548,1551 decreases. Tfiis effect is clearer wfien 
tfie same comparison is made for tfie fiigfi orbital inclination case (lower panel in Figure 
ED. We conclude tfiat for fiigfi ionization lines Civ AA1548,1551 and Nv AA1238,1242 tfie 
stronger tfie accretion rate, tfie stronger tfie line intensity. For lower ionization energy 
lines (e.g.. Si IV AA1393,1402 and Cm A1175) we found an increase in flux for lower mass 
accretion rates. 

Tfie effect of tfie primary mass is analyzed by comparing models "c" and . Tfiese 
disk models fiave, approximately, tfie same maximum disc temperature (Table [1]). However, 
because of dependence of wind parameters on Mwd? tfie system witfi lower primary mass 
(0.6 M0) fias a larger mass-loss rate (~ 6 times) tfian tfie system witfi Mwd =1 M0. At 
a fixed temperature tfie WD mass influences botfi tfie disk gravity and size of tfie disk. 
Increasing tfie WD mass increases tfie gravity at a given (r/R-wo) disc location — directly 
because of tfie increase in Mwd? and indirectly because Rwd has decreased. An increase in 
gravity translates into an increase in escape velocity, and hence a reduction in the mass-loss 
rate per unit area from a given location (as prescribed by r/Rwo) (see expression for rho, 
equation [6]). In addition, an increase in the WD mass decreases the size of the emitting 
wind region because of the dependence of T^jf on (r/Rwo) (equation [7j), which also leads 
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to a reduction in the total mass loss from the disc. Similarly, the CAK model predicts that 
the wind terminal velocity will increase with Mwd, as seen in Tabled] 

The effect is evident in Figure [HI where the ratio of mean specific intensities in the line 
absorption (I/) component over the continuum (Ic), defined as 



Ir 



is shown in the impact parameter space. In Equation fllip . the mean specific intensities 
are calculated between two frequencies that embrace the blueshifted absorption profile {voo 
z/jj, and two frequencies that embrace the nearby continuum (z/q^, Vf^). Note that the 
absorption region (Ii/Ic<l) is larger when the primary mass is lower, and consequently, the 
absorption profile is deeper as it is shown in the right panel. The same behavior is found 
for the other UV resonance lines. The intensity of the emission component depends mainly 
on the ionization state. For high ionization species, like Civ and Nv, the intensity does 
not correlate with the primary mass. For the low ionization species and Siiv AA1393,1402 
the emission component appears stronger by increasing Mwd- These lines are formed in the 
cooler external wind regions that are strongly influenced by the outermost wind radius rj 
and the external aperture angle ^2- 

We introduce a direct change in the local value aiming to study the effects that it 
could produce on the line profiles. This could be used to understand a more complex wind 
structure. The physical and geometrical parameters of model "/;." were used with a wind 
mass-loss rate increased five times (to Mit,=2x 10^^*^ M0 yr^^) for the comparison model. 
Figure [9] shows the effects of that density enhancement on line profiles for two inclinations. 
The main effect is a rise in the line emission. For i=30° there is also a slight increase in the 
depth of the absorption profiles — the profiles shown for the i=70° do not show absorption 
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components.. These effects are expected from the increased wind density. The temperature 
and ionization structures are not much affected by the enhanced mass loss rate. Those hues 
formed in the wind from inner disk regions are more affected by variations in as these 
regions have stronger winds. In our models the extension of the line wings is unaffected by 
the mass loss. 

The effect of the orbital inclination on the synthetic spectra was extensively explored. 
In Figure [10] the synthetic spectra are shown for model "d" with several inclinations. The 
inclination angle goes from 10° to 80°. The overall behavior follows the observational data 
for CVs in the UV. Low inclination models present line profiles that are preferentially in 
absorption. The atmospheres are seen almost face on, the profiles are blueshifted due to 
the strong vertical acceleration, and the line profiles are narrow because the influence of 
the wind rotation is weak. In the case of intermediate inclinations, emission structures 
begin to emerge in the profiles because high temperature wind regions are now frontwards 
of the cooler disk continuum. These emission structures produce P Cygni profiles in several 
lines such as Cm A1175, Siiv AA1393,1402 and Nil A1085. The emission structures appear 
blueshifted in intermediate inclination models because the wind regions with positive 
projected velocities are occulted by the disk, especially in the case of high ionization lines 
that are produced in the inner disk regions. When the orbital inclination increases, the 
continuum source aspect decreases relative to the wind dominated region. In this sense, for 
high orbital inclination models the continuum dims as the emission line intensity grows as 
shown in Figure [TOl The line width is also increasingly dominated by disk- wind rotation in 
high inclination models, where the emission lines become broader and the blending with 
photospheric profiles is stronger. The U V continuum inten sities strongly depend on disk 



inclination via the limb darkening effect (jPiaz et al 



19961 ). When the mass-loss rate is 



strongly increased a higher reddening is found, but this effect is weaker in the case of high 
orbital inclinations. 
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All models show emission in Si ill A1206 even in low inclination disks. Also, emission 
components were found in low inclination models for Si IV AA1393,1402. These emission 
components are not found in most low inclination CVs. In our models these features 
are mostly produced in the outermost wind regions, from gas that is coming from cooler 
portions of the disk (r^rj). The wind electron temperature there is ~11700 K, and the 
dominant ion of silicon is Si ill but with a fraction of Si iv. The outer disk atmospheres 
are larger in these models, and that is the reason why the emission is present even for low 
inclinations. A similar situation happens for Cm A1175, but in this case the Cm fraction 
in the outermost atmosphere is lower and the corresponding emission structures are weaker. 

The general behavior of the observed line profiles with the orbital inclination is well 
reproduced by models. However, these models yield some profile structures that are not 
commonly found in the UV lines of CVs, like strong Siiv AA1393,1402 emissions and a 
blueshifted emission for Civ AA1548,1551 and Nv AA1238,1242 in low inclination systems. 
These features were found to be dependent on the wind geometry and ionization structure 
that bears different emission volumes for each ionization state. 

Aiming to study the influence of an extended low temperature wind filling part of the 
primary Roche lobe, extra models were calculated with the atmosphere extending beyond 
the outermost region of the original model. This additional atmosphere wind model has a 
lower temperature and its structure is also extended until Kdisk and limited by 62 in the way 
explained in ^|2j The presence of this atmosphere shows a significant effect on the synthetic 
spectrum. The major changes occur in the low ionization lines (Cm A1175, Sim A1206 
and Cll A1335) and also in Silv AA1393,1402. These lines have their emission structures 
weakened while their absorption profiles become deeper. The high ionization lines (eg. C iv 
AA1548,1551) are not affected by the inclusion of the cooler wind models. This simulation 
reveals the importance of the outer wind ionization structure, and helps us to constrain the 
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parameters that characterize the wind and disk. 



4. COMPARISON WITH OBSERVATIONS 

In this section we confront our synthesized spectra with HST UV observations of the 
Nova-hke CVs RW Tri and the lUE data for V347 Pup. Through this comparison, the 
capabihties and hmitations of our method are revealed. The need for some model changes 
became evident from this comparison, and consequently some additional calculations were 
made. The new models were labeled as "6"' , "/i"' , "6"" and "/i"" . The model parameters are 
shown in table [6l 



4.1. RW Tri 



The Nova-like RW Tri is an eclipsing system. The orbital period is ~0.25 days 



(~5.5 hours), the pr i mary mass is arou nd ~0.55 Mq and the orbital inc 



Bitter fc Kolb 



2003 



Poole et al. 



20031 ). The distance was estimated by 



ination is ~ 70° 



VlcArthur et 



1999 ) to be between 310 and 380 pc. From a UV continuum analysis. 



al. 



Puebla et al. 



(120071 ) calculated a value for the mass accretion rate of ~ 4.6x 



is compatible with those calculated by 



Groot et al. 



(120041 ) and 



0-^ yr-^ T 



lis value 



Mizusawa et al. 



torn . A 



summary of binary system parameters is shown in Table O The UV as well as the optical 
spectrum show strong emission li nes. Observations in the UV also show that the lines are 



less eclipsed than the continuum (ICordova fc Mason 



1985 



Drew fc Verbunt 



19851). 



^T0) 



The spectroscopic data was obtained from the HST data archive (MASTj^). All spectra 
were taken by the GHRS spectrograph in RAPID mode, using the grid G140L with a 



^http:/ /archive, stsci.edu 
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spectral coverage betyeen 1 150 A and 1660 A in low resolution (~1 A) (for more details of 



data see 



Mason et al 



Cardelli et al. 



19971 ). The data were corrected fo r interstellar extincti on using the 



( 119891 ) law and a value of E{B - V)=0.1 f lBruch fc Engel 



19941 ) ■ 



In this work we used the basic parameters of the "6" and "/i" models from Table [T] in 
order to fit the RW Tri UV data. The "/i" model has physical parameters close to those 
previously estimated for RW Tri. Unlike, model "6" has a white dwarf mass higher than 
reported for RW Tri. We have also analyzed this case because this high mass model is 
capable to produce stronger C IV and N V emission. It was necessary to change some other 
parameters in order to fit the emission lines. In particular, the geometrical parameter 62 
was adjusted in order to match the observed intensity of Cm A1175 and Siiv AA1393,1402 
lines. The 9i value is used for matching the Nv AA1238,1242 and Civ AA1548,1551 lines. 
Also, in order to achieve better line ratios, we have arbitrarily included denser wind regions 
as explained in section 13.2. H by enhancing the local mass loss rate five times between 
certain radii (r^^ and r^j). 

In the top panel of Figure ITT] the UV data for RW Tri is shown together with the "best" 
calculated synthetic spectra using the "6" model and the modified model "6"'. Synthetic 
spectra were computed assuming an orbital inclination of i=70° and inner aperture angle 
61=5° and 6*2=35°. The outermost radius for the wind atmospheric structure is rf=17 Rwd 
(~9. 05x10^ cm) with Te//(r/)~16800 K. Both models in Figure [TTb show the synthetic 
spectra for a model with rc=l. 87x10^ cm. The models were scaled to a distance of 550 pc. 
We set rc lower than rj because if rc is increased to rj (~9. 05x10^ cm) the continuum flux 
is slightly increased, damaging the line-continuum ratio. The derived distance seems too 
high when compared with the distance estimates found in the literature (~350 pc). This is 
possibly due to the Mwd used in the models which is too high considering the range found 
in the literature. The modified model "6"' has a density enhancement between 1.56 Rwd 
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and 2.32 Rwd (0.83-1.25 xlO^ cm) set by increasing m(r) to 5 times the value given by our 
wind acceleration prescription. The figure shows the improvement on line emission ratios 
from model "6" (without dense region) to model "6"' (with dense region). The intensities of 
the Civ AA1548,1551 and Nv AA1238,1242 have been enhanced, as expected, but still are 
weaker than those observed. A significant improvement was found for the Si IV AA1393,1402 
line, that shows a better line ratio regarding CiV and Nv in model "6"'. 

The lower panel of Figure [11] shows the best spectra using the model "/;." and the 
modified model "/;,"'. The inner and outer aperture angles are 6i=2° and 6'2=10°. This 
model has a value closer to that estimated for RW Tri yielding a scaled distance of 

~450 pc, closer to the literature values. In the "/i"' model a dense region was again included 
to improve the lines as shown in Figure [TTl in the same sense as explained for the "6"' 
model. In the figure, the values between parentheses are the positions (in ring numbers) 
of the density enhanced region (table |6]). The models h's are cooler than b's models so 
their continua are redder, but still bluer than the observed continuum. We also found an 
improvement in the CiV AA1548,1551 and N v AA1238,1242 line intensities in h' model, but 
for the case of C iv, it was not enough to match the observations. 

In Figure [12] the emission profiles for the main UV lines are shown together with the 
better synthetic profiles (models b' and h'). We found a good agreement for the Cm A1175 
and Si IV AA1393,1402 profiles. The model with higher Mwd displays better agreement with 
the line widths. Both models show a lack of fiux in the CiV AA1548,1551 line, even when 
the inner wind regions are enhanced. This indicates that the C iv emission region is larger 
than predicted by our models. We conclude that the boundary layer emission may play an 
important role in the formation of such a strong CiV line. Model Nv AA1238,1242 lines are 
wider than the observed profiles showing that the N v line in the system is possibly coming 
from a region with a lower rotation velocity. The model "6"' line intensity is lower, as in 
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the case of C IV, but it is stronger in model h' which has a more extended (in radius) dense 
region (see table [H])- It is clear from this that the high ionization region in the wind should 
be larger. The models with M-wd =0.6 Mq show a narrow blueshifted emission structure for 
Cm A1175 and Siiv AA1393,1402 lines. This component comes from the innermost regions 
with high vertical acceleration. 

Table [7] shows the equivalent widths (EWs) for model "6"' , and data. An 
evaluation of the model versus observed profiles is also given, on the basis of the integrated 
residuals and visual inspection. We can see that model "6"' shows better values of EWs 
and also a better profile matching for Cm A1175 and Silv AA1393,1402. On the other 
hand, the high ionization lines Nv AA1238,1242 and Civ AA1548,1551 are not so good. 



Recently, 



Noe 



jauer et al 



version of the 



(120101) mo deled the UV lines of RW Tri using an improved 



metal 



Lucy 



20oJ 



(120051). Their models 



20031 ). They found 



Long fc Kniggd (l2002l ) code developed by 
incorporate the formation of recombination lines (see also 
good agreement for C iv and Si iv line intensities, but with a mass-loss rate almost 5 times 
larger than ours. Although this difference their and our values are compatible with the 
hydrodynamical models predictions. Their 3D Monte-Carlo model has a reduced Silv 
emission volume, yielding a better match of the Si iv and C iv flux ratio. They also obtain 
a strong C IV emission without include dense regions. In their models C iv fills an extended 
region in the wind, N V is limited to inner wind regions. Si IV and C ill and lower ionization 
species are limited to lower wind regions (close to disk surface) at larger radii (see fig 8 in 



Noebauer et al. 



mm)- 



In figure [13] we present the ion density structure calculated with using the model "6"' 
parameters. The strong dependence of the ionization structure on the temperature of 
radiation that is coming from disk is evident. In our approximate geometry this radiation 
arises only from disk regions that are strictly below of each wind point. That is the reason 
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why the ionization structure shows larger variations in the radial direction than those 



found in t he vertical profit 



Distinctly, 



Noebauer et al 



es. In our model structure there is no radial radiation transfer. 



(120 lOf ) with their Monte Carlo method obtained a 3D treatment 
of the extended wind regions assuming a flat disk surface. These authors found a more 
extended C IV emitting region thanks to the irradiation of the distant wind by the inner 
disk. 



4.2. V347 Pup 



The luminous CV V347 Pup is a high orbital inclination (i~80° iRitter &: Kolbl 120031 ) 
Nova-like system, with bright optical and UV emission lines. Its orbital period is ~0.23 
day s (~5.5 hours), the pr i mary mass is ^0.63 Mfn an d the distance i s estim ated as ~510 



pc (IDiaz fc Hubeny 



1999 



Thoroughgood et al 



20051) 



Puebla et al 



(120071 ) estimated a 



mass accretion r ate of ~6xl0~^ Mq yr~^, a value compatible with that calculated later by 
Ballouz fc SionI (120091 ). A summary of the binary parameters is shown in Table [5l The UV 
data were obtained from the lUE archive and were collected with the SWP camera in large 



aperture mode with a spectral resolution of ~6 A. The spectra extend from 1 



A. The data were corrected for interstellar extinctio n using the 



and a value of E{B — V)=0.06 (IMauche et al. 



Cardelli et al. 



50 A to 1950 



19941 ) 



( 119891 ) law 



In this case we used the modified models "6"" and (see table EJ- Figure shows 
the comparison of model spectra with the V347 Pup lUE data. All models have the same 
geometrical parameters: i=80°, 6i=5°, 02=4:5° (including b and h models). As is the case 
for RW Tri, a dense region was included in order to improve the intensities of the high 
ionization lines. Following the previous models, this region is five times denser than the 
original model and is located between ring 2 and 3 (0.83-1.25x10"^ cm when Mwd =1 
and 1.25-1.73x10"^ cm when M-wd =0.6 Mq). However, for this case we did not find any 
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relevant differences with the non-altered models [b and h). We only found a slight increment 
in Nv AA1238,1242 and Civ AA1548,1551 lines. The upper panel model in Figure [Til was 
scaled to 710 pc to match the flux level data in the 1450 A continuum region. This model 
shows stronger lines than those observed for Cm A1175 and Siiv AA1393,1402 lines, and 
weaker lines for Nv AA1238,1242 and CiV AA1548,1551. When using a primary mass of 
0.6 Mq (lower panel in the figure) the model intensities of Silv and Cm lines improve, 
while the intensities of N v and C iv lines decrease. This synthetic spectrum was scaled to a 
distance of 570 pc, a distance closer to that found in the literature, due to the lower Mwd- 

Figure [15] shows the line profiles of the models along with the data. A fair agreement 
for the Cm A1175 and Si IV AA1393,1402 profiles is found, but the models underestimate 
the flux in the high ionization lines (Nv AA1238,1242 and Civ AA1548,1551). This model 
behavior can be explained by the temperature structure and by the emission volume. 
The density increment at the inner wind zones helps to enhance the high ionization lines 
emission, but not enough to match observations because of the limited volume where the 
high temperature radiation ionizes the ejected gas. The enhanced intensity in C ill and 
Si IV in the model with high Mwd is due to the larger emission volume of the outermost 
wind atmosphere, where these lines are mainly produced. 

In Table M we compare EWs for models "fe"" and "/?,"" with observation. Similar to 
the RW Tri case. Cm A1175 and Si IV AA1393,1402 are well reproduced by the models 
while Nv AA1238,1242 and Civ AA1548,1551 are not. Again, we found a lack of flux for 



the last lines due to their small emission volume. [Shlosman et al.l ( Il996[ ) modeled the C iv 
AA1548,1551 line profile in and out of eclipse. They could fit the line profiles, but needed a 
very high mass-loss rate of ~lxlO~^ M0 yr~^, almost half of the accretion rate, and used 
system para meters (Mwn and Mn) that produce a hotter disk. We modeled a cooler disk 



than that of 



Shlosman et al 



( Il996[ ) (as they suggested) but could not attain the data flux 
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line. It is possible that the inclusion of other ionization sources, such as the BL emission, 
could improve the match to the Nv AA1238,1242 and CiV AA1548,1551 emission profiles. 

In figure [16] we present the ion density structure calculated using the model "/i"" 
parameters for Nv, Civ, Cm and Siiv. In this case the ^^cylindrical" structure is also 
evident. In this cooler disk, the regions with dominant Nv and Civ are closer to disk 
axis than thev are in high white dwarf mass case case, resulting in weaker lines, unlike the 



Noebauer et al. 



(120 lOl ) models where, specially for C IV, these regions are fairly extended 



through the wind. 



5. DISCUSSION 

In our simulations, the vertical wind velocity is calculated taking into account the local 
physical environment for each wind atmosphere. Previously, the disk and wind were treated 
independently or with a fiducial transition region. Also, in earlier works, the ionization 
equilibrium and level populations have been calculated using simplified methods. Even 
with several simplifying assumptions, we found that our model spectra resulting from our 
complete synthesis present emission lines and continuum fluxes close to those observed in 
bright CVs. By comparing observed data of two high inclination systems with the models 
we have found similarities between the model and actual structured line profiles. However, 
we include denser regions within the wind in order to obtain reasonable line intensities, 
in particular for matching the high ionization line fluxes. The presence of c ondensations 



was previously predicted by hydrodyn amical models of accretion disk winds (jPereyra et al. 



1997 



2000 



Pereyra Sz Kallman 



20031 ). The figures fTTl and show that it was not enough 
to exactly match the observed profiles. We found good agreement with observations for 
Si IV and Cm lines. For high ionization lines Civ AA1548,1551 and Nv AA1238,1242 the 
model lines are weaker than those observed. This could be due to the infiuence of outer 
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wind irradiation by the inner disk and/or from the effect of the boundary layer on the 
ionization structure, which are not yet considered in the models. 

By analyzing the models for low orbital inclination systems, we found that the synthetic 
line profiles do not adequately reproduce the line profiles. However, for such cases we found 
that the profiles are sensitive to the wind geometry and vertical structure. The vertical 
structure is affected by the gas expansion and geometric radiative dilution that are not 
taken into account in the plane-parallel approximation. In addition, the interaction among 
neighboring atmospheres must be taken into account. This may be crucial for CVs — here 
we tried to solve a 3D (or at least 2.5D) problem using a set of ID models. Nevertheless, 
the models reproduce to reasonable accuracy the high orbital inclination system data. This 
leads us to the conclusion that the line profiles of high inclination systems are less dependent 
on the geometry of the wind, being dictated mostly by the accretion disk parameters. 

In previous works the attention often focussed on CiV AA1548,1551 , because of its 
strong sens itivity to wi nd parameters. Nevertheless, it was shown (e.g. 



fl2002h and 



Long fc Knigge 



Longl (120061 )) that is an arduous problem to reproduce other UV lines with the 



same wind parameters. In the present models we attempt to fit simultaneously several line 
intensities and their profiles. Our results suggest that a highly structured wind can help to 
improve the fit to the line profiles and the intensity ratios. In addition, we found that it is 
important to have a 2D consistent structure that take into account the influence of disk 
radiation on the whole wind. We found model absorption pro files similar to those observed 



in the FU SE spectral region (ILong 



2006 



Froning et al. 



20031 ). but with blue shift velocities 



higher than observed. This may suggest that these lines come from the photosphere wind 
transition region, but the wind is less accelerated than calculated here. 



In a recent wor' 



by 



c. iNoebauer et al.l ( 120101 ) uses a 3D Monte Carlo method developed 



Long fc Kniggd ( 120021 ). They show more complex ionization structures, where Civ 
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dominates large extensions of the wind, which bears a stronger emission in the CiV 
AA1548,1551 hne. We also found an extended CiV region in our winds, but highly 
concentrated close to the disk axis. 

On the other hand, our method is capable of take account consistently the photosphere- 
wind transition region, although for that we had to give up the 2.5 D description for the 
wind temperature. We have started a brief analysis of the effects of that region on the 
emitted spectrum. A diagnostic feature of such regions are the blueshifted absorption 
components found in low orbital inclination models. We tried to analyze the effect of that 



interface region on the emitted spectrum using the Cll A1335 



ine. T his is a photospheric 



20101). We used two models 



line or it is produced in the wind close to disk (iNoebauer et al.l 
with different acceleration laws and an orbital inclination 2=30° (models h and i from table 
[1]) in order to test the effect on the profiles. The figure [17] shows that even considering 
that the line is produced very close to disk photosphere, the acceleration law has a strong 
influence on the profile and on its equivalent width. These facts show that the interface 
region structure is important to analyze the spectrum. 

Another obstacle in the comparison of disk wind models with observations is the 
variability seen in the UV resonance lines. Strong varia b ility, n ot related to the orbital 



phase, was already found in CV wind lines. 



Prinja et al 



(l2000bl ) showed that the nova-like 



BZ Cam has rapid line variability, on timescales as short as minutes suggesting the existence 
of perturbations that could travel through the wind on th at time scale. Other systems show 



variability in differ ent time scales, for in stance V603 Aql (iPrinia et al. 



( jPrinja et al. 



20031). On the other hand. 



Hartley et al. 



2mm and RW Sex 



(120021 ) found variability patterns in 



the high accretion rate systems IX Vel and V3885 Sgr that could challenge the line-driven 
hypothesis for those winds. 



The wind emission has been pointed out as a candidate for solution to the known 
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the c olor-magnitude problem seen in previous models (e.g. IWadd Il988l : iKnigge et al. 



19981 ). We found a slight reddening of the spectrum due to the wind which is stronger 
for higher mass-loss rates and A < 1500 A. This is due to the height where the continuum 
becomes transparent, which is slightly higher when the mass- loss rate is increased, and 
the corresponding temperature is lower. Howev er, this redde n ing is not sufficient to fully 



overcome the UV continuum emission problem. iPuebla et al.l (120071 ) found an average color 
mismatch of ~0.3 mag between 1500-3000 A. Our dense model lowers that value to ~0.24 
mag when the 1000-3000 A region is taken into account. 

It is important to summarize the main limitations and the scope of our approximations. 
The strongest limitation is not taking into account the radiation of the whole disk on the 
local wind structure (temperature, electron density, population levels and radiation itself). 
This is more significant for the lines that are formed high above the disk. This point will 
be improved in future work by adding one dimension to the structure calculus. 

The outer (cooler) atmosphere is mostly responsible for the low ionization emission 
found in our models that are not observed in low inclination systems. Again, in our 
approximation a possible explication is that the outer and higher wind regions do not see 
the radiation that is coming from the hotter inner disk (and boundary layer as well). Hence, 
a large low ionization volume dominates the outer regions of the wind. 

We arbitrarily included regions with increased wind density trying to improve the line 
fitting. This also explore the possibility of more complex structures that could produce 
the line emission. This is not intended to be a definitive or unique model and further 
exploration is needed. 

However, our synthetic spectra present several trends that can be contrasted with 
observations, for example: i) the dependence of the high-to-low ionization line ratios with 
the wind temperature, ii) the dependence of the depth of the absorption components with 
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the WD mass, iii) the emission hne strength as a function of the disk temperature, iv) our 
models are able to reproduce the observed rise in the emission lines with orbital inclination 
and increasing mass-loss rates. The UV FOS and lUE data for the high inclination systems 
LX Ser, V348 Pup and V347 Pup (systems with different disk temperature) show some 
of these trends: hotter disks produce stronger lines, and the line ratio Civ/SilV falls for 
cooler disks. 

As was already pointed out, a new ingredient of the models is the photosphere wind 
interface. The effect of such a region on the synthetic spectrum is evident mainly in the 
photospheric lines (e.g. Cll A1335) and in the blushifted absorption features. In the later 
case, however, this effect can be overshadowed by the absorption in the upper regions of the 
wind. The photosphere-wind interface is insensitive to the of radiation produced at distant 
regions in the disk. 

The next step to improve our models is to abandon the cylindrical geometry generated 
by the ID atmosphere structures, trying to use a hydrodynamically computed geometry to 
reduce the numbers of free parameters. This alone may not solve all the problems, but will 
facilitate the model constraining itself. In addition, a method to include the influence of 
the whole disk radiation on the wind structure will be developed, possibly with the help of 
Monte-Carlo techniques. 

6. SUMMARY AND CONCLUSIONS 

In this work we made an attempt to understand the UV accretion disk and its 
wind emission. We use a simple 2.5D model formed by a set of ID models that treat 
disk-wind interface self-consistently. In general, the geometry of CV systems does not 
allow a one-dimensional treatment of the problem. However, in order to make the problem 
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tractable, we constructed a grid of ID models based on local disk properties, and then 
combined these models into a 2.5D model. Disk rotation was taken into account — this is 
essential for computing line profiles. 

The density and velocity structure of the wind above the disk was determined 
using line-driven wind theory. Due to the dependence of gravity on height, the wind 
acceleration laws differ from stellar results. Both the atmospheric structure and radiative 
transfer depend on the wind acceleration. We found that the ionization structure and 
wind temperature strongly depend on the vertical velocity. The Euler equation solutions 
yield strongly accelerated winds, with maximum acceleration close to the photosphere. 
Because of this, the atmospheric structure shows strong variations in the photosphere-wind 
transition region. Higher in the wind, the physical structure is almost constant with z, and 
the ionization state is close to that corresponding to a stellar atmosphere with temperature 



From an analysis of synthetic spectra, we find that the wind spectral features are 
strongly affected by the classical disk properties. The disk temperature sets the wind 
ionization structure in the sense of the size and location of line formation regions. The 
primary mass affects the strength of absorption profiles as well as the line emission 
volume. An increase in the mass loss rate (and wind density), as expected, strengthens 
the wind characteristics of profiles; the absorption components are deeper and the emission 
structures are enhanced. An increase in the mass-loss rate also makes the continuum 
slightly redder especially for A<1500 A. However, this reddening is not enough to overcome 
the color-magnitude problem which remains present in our models. In agreement with 
observations, absorption and P Cygni lines are synthesized for low inclination disks while 
an emission line spectrum is found for high inclination disks. 

While many characteristics of CV spectra can be explained by our simplified models. 
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discrepancies remain. In future modeling it will be important, for example, to take into 
account the interaction between different regions of the disk wind. In order to reproduce 
the UV line intensity ratios observed in high inclination CVs, and to explain variability, 
condensations in the wind may be necessary. 

We thank FAPESP (process: 2005/04128-5) for financial support. MPD also 
acknowledges support from CNPq under grant ^^305725. DJH acknowledges partial support 
by the National Aeronautics and Space Administration under grant ATP 03-0104-0144. We 
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Fig. 1. — Geometry of system used to model the disk wind and to calculate the spectrum 
synthesis. The geometrical parameters for spectral synthesis are thus defined (see text). 




0.0 0.2 0.4 0.6 0.8 1.0 



Fig. 2. — Combined and velocity fields. Red lines show the radial assembly of wind 
atmospheres. Blue line shows a hyperbolic stream line with an asymptotic aperture angle /3. 
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0.5 LO 1.5 2.0 

r [10^° cm] 

Fig. 3. — Euler equation solutions for a disk wind with Mwd=1 and Ma=10~^ M© 
yr~^. The upper panel shows the temperature distribution. The middle panel shows the 
terminal velocities for the upper branch solutions for the three a (0.9, 0.7 and 0.5) values 
compared with escape velocity. The bottom panel shows mass loss fluxes that result from 
critical solutions. 
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r [10^ cm] 

Fig. 4. — Radial mass loss contribution in Mq yr~^ (black line) and mass loss flux distribution 
in gr s~^ cm~^ (red line). Both correspond to model "e" (table [3]). 




Fig. 5. — Vertical temperature profiles for wind atmospheres. All rings in model "e" are 
shown. 
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Fig. 6. — Carbon vertical ionization structure for ring 6 of model "e" (Te//=30300 K). The 
left panel shows the inner atmosphere, photosphere-wind transition. The right panel shows 
the extended wind region. 
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Fig. 7. — Effect of accretion rate on the line profiles for Cm A1175 and CiV AA1548,1551. 
Ma = l.OxlO-s Mq yr-i (solid line), S.OxlQ-^ M© yr"^ (dotted line), l.OxlQ-^ M© yr"^ 
(dashed line) and 5.0 x 10^^° M© yr"^ (dot-dashed line). The orbital inclination for all models 
in the upper panels is i=30° and 2=70° in the lower panels. 
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P (10" cm) velocity [km s '] 

Fig. 8. — The ratio of the mean specific intensities in the hne and in the continuum in 
the spectral region of CiV AA1548,1551 for models "c" and "/i" is shown for an observer's 
inclination of 30°. Only the absorption regions {h/Ic < 1) are shown. The right panel shows 
the calculated line profiles obtained by integrating across the whole disk. The models have 
similar disk temperatures but different primary masses, 1 (solid line) and 0.6 (dotted 
line). 
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Fig. 9. — Effects of the wind mass-loss rate on tlie line profiles of Siiv AA1393,1402 and 
Civ AA1548,1551. Synthetic profiles are shown for model "/i" with M^=4 (solid line) and 
M^=20 (dotted line), both in 1x10^^^ Mq yr~^ units. Both models were calculated with an 
orbital inclination of i=30° (upper panel) and i=70° (lower panel). 
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Fig. 10. — Mid and far ultraviolet disk model spectra at several orbital inclination angles for 
model "(i". The angles span the values from 10° to 80°, with 10 degree steps. The spectra 
are scaled to a distance of 400 pc and separated by a constant offset of 1.5x10"^^ erg cm~^ 
s— 1 for clarity. 
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Fig. 11. — Best simulations for RW Tri UV data. The data are the mean of spectra taken out 
of echpse. The upper panel shows the best-fit models using the disk parameters for model 
"fe" and "6"' (see text for details). The lower panel shows the best-fit models using the disk 
parameters for model "/;." and "/;,"'. The numbers in the parentheses show the rings where 
a density enhancement was included for the altered models (see text for details). 
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Fig. 12. — Continuum normalized line profiles for the main UV lines of RW Tri: Cm A1175, 
Nv AA1238,1242, Siiv AA1393,1402 and Civ AA1548,1551. The parameters for each model 
are given in the text. 
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Fig. 14. — Best simulations for V347 Pup UV data. The data are the mean of a spectra taken 
out of echpse. The upper panel shows the best model using the disk parameters for model 
b and b". The lower panel shows the best-fit models using the disk parameters for model h 
and h". The numbers in the parentheses show the rings where a density enhancement was 
included (see text for details). 
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Fig. 15. — Continuum normalized line profiles for the main UV lines of V347 Pup: Cm 
A1175, N V AA1238,1242, Si IV AA1393,1402 and Civ AA1548,1551. The parameters for each 
model are given in the text. 
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Fig. 17. — C II A1335 line profiles for models h and i and an orbital inclination of 30°. They 
show the effect of the different acceleration laws on the photosphere-wind interface region. 
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Table 1. Model Parameters 



Model 


MwD 


RwD 


Ma 




M„ 




T 1 

disk 


V 2 


V 3 




a 


k 






Mq 


10^ cm 


M© yr " 


-1 


Mq yr - 


-1 


10* K 


km s^^ 


km s^^ 


10^ cm 








a 


1.0 


5.5 


1.0x10" 


-8 


1.1x10" 


-10 


8.60 


1540^ 


491 


13.9 


0.9 


0.6 


Ci 


b 


1.0 


5.5 


5.0x10" 


-9 


4.0x10" 


-11 


7.20 


1640^ 


461 


9.0 


0.9 


0.6 


Ci 


c 


1.0 


5.5 


1.0x10" 


-9 


7.0x10- 


-12 


4.80 


1470^ 


580 


6.2 


0.9 


0.6 


Ci 


d 


1.0 


5.5 


5.0x10- 


-10 


5.0x10- 


-12 


4.10 


1110^ 


570 


4.2 


0.9 


0.6 


Ci 


e 


0.8 


6.7 


1.0x10" 


-8 


9.3x10" 


-11 


6.80 


6200^ 


1200 


13.9 


0.7 


0.7 


C2 


f 


0.8 


6.7 


1.0x10" 


-8 


4.7x10- 


-10 


6.80 


5800^ 


908 


13.9 


0.5 


0.9 


C2 


9 


0.8 


6.7 


1.0x10" 


-9 


1.7x10"" 


-12 


3.86 


1248^ 


324 


4.9 


0.9 


0.6 


Ci 


h 


0.6 


8.6 


5.0x10" 


-9 


4.0x10"" 


11 


4.63 


1182^ 


523 


9.7 


0.9 


0.6 


Ci 


i 


0.6 


8.6 


5.0x10" 


-9 


6.5x10- 


11 


4.63 


4500^ 


1360 


9.7 


0.5 


0.9 


C2 



^Maximum disk effective temperature (Eq. [7]). 

^Maximum vertical velocity V^^ (Eq. [S]) for the first ring. 

■^MaximmTL vertical velocity Vz^ calculated with Eq. [5] for the last ring at r/. 

^External wind atmosphere radius, T(rj^)^14000 K. 

^Lower branch solutions. 

^Upper branch solutions. 
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Table 2. Geometrical Parameters 

Parameter Values 

i (°) 30, 40, 50, 60, 70, 80 

ri/RwD 1-05 

Oi n 0, 5, 15 

62 (°) 0, 15, 45 

Zmax/Rrfisfc I'O 

rc/rf 1.0 
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Table 3. Physical characteristics of each ring used to calculate the model "e", namely: the 
radius of each ring, the mass loss flux, the terminal velocity and effective temperature. The 

star-like wind parameters are also listed. 



Disk parameters Star-like parameters 



Ring 


r 

RwD 


rri(r) 

gr cm^^ s 


-1 


Woo 

km s"^ 


Te// 

10"* K 


L 

L© 


M„ 
Mq yr-i 


1 


1.05 


2.6x10" 


-5 


4043 


5.42 


7.9x10'"^ 


1.2x10"*^ 


2 


1.52 


1.4x10" 


4 


3351 


6.87 


2.0x10^ 


5.5x10"'^ 


3 


2.20 


9.6x10" 


■5 


2770 


5.96 


1.1x10^ 


3.6x10"^ 


4 


3.19 


4.6x10" 


5 


2500 


4.86 


5.1x10^ 


1.7x10"^ 


5 


4.63 


1.9x10" 


5 


2500 


3.86 


2.0x10^ 


8.3x10"'^ 


6 


6.70 


7.3x10" 


■6 


2400 


3.03 


7.7x10"* 


3.2x10"'^ 


7 


9.72 


2.7x10" 


■6 


1950 


2.35 


2.7x10"* 


1.2x10"'^ 


8 


14.08 


9.6x10" 


■7 


1600 


1.81 


9.9x10^ 


9.3x10-8 


9 


20.41 


3.3x10" 


■7 


1260 


1.39 


3.4 xlO^ 


1.4x10"^ 


10 


29.58 


1.2x10" 


■7 


980 


1.06 


1.19x10^ 


1.1x10"^ 


11 


42.80 








0.81 






12 


62.10 








0.62 
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Table 4. Ionization states used to calculate the atmospheric structures. 



Ion 


Levels 


Super-levels 


Ion 


Levels 


Super-Levels 


H I 


20 


20 


V 


41 


41 


Ho T 
lie i 






yj V i 




1 '? 


He II 


22 


22 


VII 


1 


1 


C II 





14 


Ne II 





14 


C III 


62 


62 


Ne III 


23 


13 


CIV 


13 


13 


Ne IV 


17 


17 


N II 


36 


36 


Si III 


20 


20 


N III 


34 


34 


Si IV 


22 


22 


N IV 


90 


90 


Fe III 


135 


135 


N V 


33 


33 


Fe ml 


96 


96 


N VI 


1 


1 


Fe IV 


100 


100 


II 


46 


46 


Fe V 


139 


139 


III 


92 


92 


Fe VI 


44 


44 


IV 


29 


29 


Fe VII 


29 


29 



1 Ions used when T^ff < 20000 K. 



Table 5. System Parameters 



System Porb Mwd M2 d E{B-V) References 

days M0 M0 pc Mq yi"^ 

RWTri 0.25 0.55 0.35 310-380 4.6x10-^ 0.1 1,2,3,4 

V347 Pup 0.23 0.63 0.55 510 6.0x10^^ 0.06 1,4,5,6 



19991 . (4) 



Refere nces. — (1) iRitter fc Kolbll2003l. (2) IPoole et al. 



Puebla et al. 



20071 . (5) 



Thoroughgood et al. 



2005 



(6) 



20031 f3nMcArtliur et al 



Diaz &: Hubeny 



1999 
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Table 6. Models Parameters 



RW Tri V347 Pup 



Parameter 


I,' 




n 


1.11 
b 


n 


Porb (days) 


0.3 


0.3 


0.3 


0.3 


Mwd(M0) 


1.0 


0.6 


1.0 


0.6 


(Mq yr-i) 


5x10-9 


5x10-9 


5x10-9 


5x10-9 


(Mo yr-i) 


8.2x10-" 


1.1x10-1° 


8.2x10-" 


8.1x10-" 


z(°) 


70 


70 


80 


80 


Rdisk (10^ cm) 


45.9 


39.7 


45.9 


39.7 


Ti (10^ cm) 


0.58 


0.84 


0.58 


0.84 


Tf (10^ cm) 


9.05 


9.74 


9.05 


9.74 


rc (10^ cm) 


1.87 


9.74 


9.05 


9.74 


Oi n 


5 


2 


5 


5 


02 n 


35 


10 


45 


45 


- (10^ cm)^ 


0.83-1.25 


1.25-3.47 


0.83-1.25 


1.25-1.73 



Radial region of enhanced wind density. See text. 
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Table 7. Equivalent widths (EW) for two RW Tri models. 



Line 


EWjv/ 


EWb' 


EW;,, 


Profile Match 




(A) 


(A) 


(A) 


Model 


Model h'^ 


Cm 


12.3 


18.3 


15.5 


good 


good 


Nv 


25.5 


43.6 


57.0 


medium 


bad 


Si IV 


28.7 


29.0 


28.0 


good 


good 


Civ 


75.0 


73.0 


53.6 


good 


medium 



Measured equivalent width from RW Tri UV data. 
^Modified models from Table [T] standard parameters (see text). 



Table 8. Equivalent widths (EW) for two V347 Pup models. 



Line 


EWm' 


EWb" 




Profile Match 




(A) 


(A) 


(A) 


Model b"^ 


Model h"^ 


Cm 


50.0 


70.0 


57.0 


medium 


medium 


Nv 


167.0 


76.0 


62.0 


bad 


bad 


Si IV 


127.0 


137.0 


110.0 


medium 


good 


Civ 


244.0 


134.0 


121.0 


bad 


bad 



^Measured equivalent width from V347 Pup UV data. 
^Modified models from Table [1] standard parameters (see text). 



